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O-C Abstract 

' A qualitative analysis of the chiral phase transition in QCD at non-zero baryon density is per- 

formed. It is assumed that at zero baryonic density, p = Q, the temperature phase transition is of the 
second order and quark condesate r] =\ {Q\uu\ 0) | = | (0 | dd | 0) | may be taken as order parameter 
of phase transition. It is demonstrated, that the proportionality of baryon masses to quark conden- 

■ sate in the power 1/3, ms ~| (0 | gg | 0) \^^^ is valid in the wide interval of quark condensate values. 
^ ' By supposing, that such specific dependence of baryon masses on quark condensate takes place up 

, to phase transition point, it is shown, that at finite baryon density p the phase transition becomes 

OO ' of the first order at the temperature T — Tph{p) for p > 0. At temperatures rcont(p) > T > Tph(p) 

' there is a mixed phase consisting of the quark phase (stable) and the hadron phase (unstable). At 

0^ , the temperature T = Tcont(p) the system experiences a continuous transition to the pure chirally 

' symmetric phase. 
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1 Introduction 



It is well known, that the chiral symmetry is valid in perturbative quantum chromodynamics (QCD) 
with massless quarks. It is expected also, that the chiral symmetry takes place in fuU-perturbative 
and nonperturbative QCD at high temperatures, (T > 200 MeV), if heavy quarks {c,b,t) are ignored. 
, The chiral symmetry is strongly violated, however, in hadronic matter, i.e. in QCD at T = and low 
' density. What is the order of phase transition between two phases of QCD with broken and restored 
5^ , chiral symmetry at variation of temperature and density is not completely clear now. There are different 
opinions about this subject (for a detailed review see Ref. 0112] and references therein). 

In this talk I discuss the phase transitions in QCD with two massless quarks, u and d. Many lattice 
calculations Pl |31 El El indicate, that at zero chemical potential the phase transition is of the second 
order. It will be shown below, that the account of baryon density drastically changes the situation 
and the transition becomes of the first order, and, at high density, the matter is always in the chirally 
symmetric phase. 

The masses of light u, d, s quarks which enter the QCD Lagrangian, especially the masses of u and 
d quarks, from which the usual (nonstrange) hadrons are built, are very small, m„,TOd < 10 MeV as 
compared with characteristic mass scale M ^ 1 GeV. Since in QCD the quark interaction proceeds 
through the exchange of vector gluonic field, then, if light quark masses are neglected, QCD Lagrangian 
(its light quark part) is chirally symmetric, i.e. not only vector, but also axial currents are conserved and 
the left and right chirality quark fields are conserving separately. This chiral symmetry is not realized 
in hadronic matter, in the spectrum of hadrons and their low energy interactions. Indeed, in chirally 
symmetrical theory the fermion states must be either massless or degenerate in parity. It is evident, 
that the baryons (particularly, the nucleon) do not possess such properties. This means, that the chiral 
symmetry of QCD Lagrangian is not realized on the spectrum of physical states and is spontaneously 
broken. According to Goldstone theorem spontaneous breaking of symmetry leads to appearance of 
massless particles in the spectrum of physical states - the Goldstone bosons. 
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Let us first consider a case of the zero baryonic density and suppose tliat the phase transition from 
chirality violating phase to the chirality conserving one is of the second order. The second order phase 
transition is, generally, characterized by the order parameter ?/. The order parameter is a thermal average 
of some operator which may be chosen in various ways. The physical results are independent on the choice 
of the order parameter. In QCD the quark condensate, rj = |(0|mu|0}| = |(0|(i(i|0}| > 0, may be taken 
as such parameter. In the phase of broken chiral symmetry (hadronic phase) the quark condensate is 
non-zero and has the normal hadronic scale, {0 \ qq \ 0) — {250MeV)^ , in the phase of restored chiral 
symmetry it is vanishing. 

The quark condensate has the desired properties: as it was demonstrated in the chiral effective 
theory [3 E], 7? decreases with the temperature increasing and an extrapolation of the curve 'r](T) to 
the higher temperatures indicates, that r] vanishes at T = Ti^-* « 180 MeV. Here the superscript "0" 
indicates that the critical temperature is taken at zero baryon density. The same conclusion follows from 
the lattice calculations (3) El E], where it was also found that the chiral condensate f] decreases with 
increase of the chemical potential |1U[ [TT] . 

Apply the general theory of the second order phase transitions jjl] and consider the thermodynamical 
potential <&(?7) at the temperature T near Since rj is small in this domain, ^(ji) may be expanded 

in rj: 

Hv) = ^o + Ia7]^ + ^B7^\ B>0. (1) 

For a moment we neglect possible derivative terms in the potential. 

The terms, proportional to rj and 77^ vanish for general reasons In QCD with massless quarks the 
absence of 77 and 77'^ terms can be proved for any perturbative Feynman diagrams. At small t = T — Tc^^ 
the function A{t) is linear in t: A{t) = at, a > 0. If i < the thermodynamical potential $(7;) is minimal 
at 7; 7^ 0, while at t > the chiral condensate vanishes 77 = 0. At small t the t-dependence of the 
coefficient B{t) is inessential and may be neglected. The minimum, 77, of the thermodynamical potential 
can be found from the condition, d^/drj = 0: 



r ^^at/B, t<0; . . 

''-\o, 00. 

It corresponds to the second order phase transition since the potential is quartic in 77 and - if the derivative 
terms are included in the expansion - the correlation length becomes infinite at T = Ti"' . 



2 Nucleon mass and quark condensate 

I show now, that the existence of large baryon masses and the appearance of violating chiral symmetry 
quark condensate are deeply interconnected and even more, that baryon masses arise just due to quark 
condensate. I will use the QCD sum rule method invented by Shifman, Vainstein and Zakharov [13], in 
its applications to baryons [14-17]. (For a review and collection of relevant original papers see [18]). The 
idea of the method is that at virtualities of order ~ 1 GeV^ the operator product expansion (OPE) 
may be used in consideration of hadronic vacuum correlators. In OPE the nonoperturbative effects reduce 
to appearance of vacuum condensates and condensates of the lowest dimension play the most important 
role. The perturbative terms are moderate and do not change the results in essential way, especially in 
the cases of chiral symmetry violation, where they can appear as corrections only. 

For definiteness consider the proton mass calculation [14,15]. Introduce the polarization operator 

n(p) = 7 J dW-(o|rc(x),c(o)|o) (3) 

where ^{x) is the quark current with proton quantum numbers and is chosen to be space-like, 

< 0, Ip^I 1 GeV^. The current f is the colourless product of three quark fields, ^ = e"''^ q'^q^q'^, q = 
u, d, the form of the current will be specialized below. The general structure of n(p) is 

^{p)^Ph{p) + h{p) (4) 

The first structure, proportional to p is conserving chirality, while the second is chirality violating. 
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Figure 1. The bare loop diagram, contributing 
to chirality conserving function /i(p^): solid lines 
correspond to quark propagators, crosses mean 
the interaction with external currents. 



Figure 2. The diagram, corresponding to chi- 
rality violating dimension 3 operator (quark con- 
densate). The dots, surrounded by circle mean 
quarks in the condensate phase. All other nota- 
tion is the same as on Fig.l. 



For each of the functions fi{p^), i = 1,2 the OPE can be written as: 

/i(f')=E ^^^H/)(0|O«|0) (5) 

n 

where (0|Oi''|0) are vacuum expectation values (v.e.v) of various operators (vacuum condensates), Cn^ 

are coefficient functions calculated in QCD. For the first, conserving chirality structure function fi{p'^) 
OPE starts from dimension zero {d = 0) unit operator. Its contribution is described by the diagram of 
Fig.l and 

Pfiiv^) ~ C'oPP'*'"[A^/(— p^)] + polynomial, (6) 

where Co is a constant, A„ is the ultraviolet cutoff. The OPE for chirality violating structure /2(p^) 
starts from d = 3 operator, and its contribution is represented by the diagram of Fig. 2: 

A2 

/2(/) = Cip'{Q\QqqWnj-^+ polynomial (7) 

Let us for a moment restrict ourselves to this first order terms of OPE and neglect higher order terms 
(as well the perturbative corrections). 

On the other hand, the polarizaion operator (3) may be expressed via the characteristics of physical 
states using the dispersion relations 

1 f Imfi{s') , „ 2 



fi(s) = — / — — ^ ds' + polynomial, s = —p^ (8) 

n J s' + s 

The proton contribution to Imll{p) is equal to 

Imn{p) = n{0\^\p)(j)\^\0)6{p^ - m^) = ttX%{p + m)6{p'^ - m^), (9) 

where 

m\p) = ^Nvip), (10) 

Xn is a constant, v(j>) is the proton spinor and m is the proton mass. Still restricting ourselves to this 
rough approximation, we may take equal the calculated in QCD expression for Il{p) (Eq.'s(6),(7)) to its 
phenomenological representation Eq.(9). The best way to get rid of unknown polynomial, is to apply to 
both sides of the equality the Borel(Laplace) transformation, defined as 

(\ n oo 
-dS) ^^'^ = n / dslmfis)e-^/^ (11) 
/ ^ 

if f{s) is given by dispersion relation (8). Notice, that 

(n-l)!(M2)"-i ^^^^ 
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Owing to the factor l/(n — 1)! in (12) tlic Borcl transformation suppresses the contributions of high 
order terms in OPE. 

Specify now tlie quark current ^(x). It is clear from (9) that proton contribution will dominate in 
some region of the Borel parameter ^ rn? only in the case when both calculated in QCD functions /i 
and /2 are of the same order. This recjuirement, together with the requirements of absence of derivatives 
and of renormcovariance fixes the form of current in unique way (for more details see [14,15]): 

^{x) = e«^^(w«C7^u^)7^75rf" (13) 

where C is the charge conjugation matrix. With the current t/(.x) (13) the calculations of the diagrams 
Fig. 2 can be easily performed, the constants Co and Ci are determined. The terms of OPE up to 
dimension 8 were accounted in the sum rule for the conserving chirality structure function fi{p^) and 
the terms up to dimension 7 were accounted in the violating chirality structure function f2{p^)- The 
perturbative corrections are neglected. The phenomenological parts of the sum rules were represented 
by contributions of proton and continuum, the latter were transferred to QCD parts. The resulting sum 
rules after Borel transformation have the form [16,17] 

M^E^{M)L-^"' + -a^L^/g + -bM^Eo{M)L-^/^- 
3 4 

1 2 2 1 T2 / '^^ > 

- 3" """M^ = ^''''^^['JP j ^^^^ 



2M4iJ,(M) + ^^4--^ 
^ ^ 81 TT 12 



wA^exp -— (15) 



M2 

a = -(27r)^(0 \qq\0) (16) 
6=(27r)2(0|^G^jO) (17) 

TT 

A" 

ff(0 I g I 0) = ml{0 I qq I 0) (18) 



= 4 / z^'e-'^dz (19) 
n! J 








Q:,,(/i) 



(20) 



Qs(M) 

= 327r2A^ (21) 

L{M, fi) in (14), (15) account the anomalous dimensions of the current ^(x) and various operators, En{j^) 
accounts the continuum contribution, which starts at s = sq. The ratio of (15) to (14) gives for proton 
mass 

m = 2a^/(M2) (22) 

The inspection of the r.h.s. of (21) at the values of QCD parameters a = 0.63 GeV'^, b = 0.24 GeV^, 
ml = 0.8 GeV^, sq = 2.3 GeV^ shows, that at 0.8 < < 1.4 GeV^ with a good accuracy /(M^) is 
equal to 1 and the value = is close to the middle in the plateau. So we get a simple formula for 
proton mass [14] 

m = [-2(27r)2(0|9g|0)]V3 (23) 

Let us now check, if eq. (23) is valid at variation of quark condensate value. Perform the scale variation 
(0 I gg I 0) ^ 7(0 I gg I 0), m ^ j^^^m, M ■y'^/^m, sq 7^^^so. Then as can be seen from the (14,15) 
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the only violations of eq.(23), come from terms, proportional to b and in (14), (15). But they are very 
small. Therefore, even at 7 = 1/3, the relation (23) changes no more, than by 7%. To go to lower values 
of 7 is not possible, because of strongly increasing perturbative corrections. 

This formula demonstrates the fundamental fact, that the appearance of the proton mass is caused 
by spontaneous violation of chiral invariance: the presence of quark condensate. (Numerically, (23) gives 
the experimental value of proton mass with an accuracy better than 10%). 

In the same way, the hyperons, isobar and some resonances masses were calculated, all in a good 
agreement with experiment [14,16]. I will not dwell on these results. The main conclusion is: the origin 
of baryon masses is in spontaneous violation of chiral invariance - the existence of quark condensate in 
QCD. 

3 The chiral phase transition in the presense of finite 
baryon density 

I would like to consider here the influence of baryon density on the chiral phase transition in hadronic 
matter. Kogan, Kovner and Tekin [19] have suggested the idea, that baryons may initiate the restoration 
of chiral symmetry, if their density is high - when roughly half of the volume is occupied by baryons. 
The physical argument in favour of this idea comes from the hypothesis (supported by calculation in 
chiral soliton model of nucleon [20]), that inside the baryon the chiral condensate has the sign opposite 
to that in vacuum. This hypothesis is not proved. Even more, it is doubtful, that the concept of quark 
condensate inside the nucleon can be formulated in a correct way in quantum theory. But the idea on the 
strong influence of baryon density on the chiral phase transition looks very attractive. For this reason 
no assumption on the driving mechanism of chiral phase transition at zero baryon density will be done 
here. The problem, under consideration is: how the phase transition changes in the presence of baryons. 
The content of my talk closely follows ref.21. 

For completeness it must be mentioned, that the variation of quark codensate with temperature is 
not the only source of baryon effective mass shift. At low T the effective baryon mass shift arises also 
due to interaction with pions in thermal bath [22,23]. However, this mass shift, which may be called 
external (unlike the internal, arising from variation of quark condensate) is related to effective mass, i.e. 
to propagation of baryon in the matter and has nothing to do with the properties of the matter as a whole 
and the phase transition. (A similar phenomenon takes place in case of vector mesons, where because of 
interaction with pions in the thermal bath the mixing with axial mesons arises [24].) 

Consider the case of the finite, but small baryon density p (by p we mean here the sum of baryon 
and anti-baryon densities). For a moment, consider only one type of baryons, i.e. the nucleon. The 
temperature of the phase transition, Tph, is, in general, dependent on the baryon density, Tph = Tph(p), 
with Tph(yO = 0) = Tc^^ At T < Tph(/9) the term, proportional to Ep, where E = y^p"^ + is the baryon 
energy, must be added to the thermodynamical potential ^ . As was shown above the nucleon mass m (as 
well as the masses of other baryons) arises due to the spontaneous violation of the chiral symmetry and is 
approximately proportional to the cubic root of the quark condensate: m = crj^^^, with c = (Stt^)^/'^ for a 
nucleon. At small temperatures T the baryon contribution to $ is strongly suppressed by the Boltzmann 
factor e~^l'^ and is negligible. Below we assume that the proportionality m ~ t]^!^ is valid in a broad 
temperature interval up to phase transition point. Arguments in favor of such an assumption are based 
on the expectation that the baryon masses vanish at T = Tph(p) and on the dimensional grounds. Near 
the phase transition point E = ^/p^ + rn? w p + rj^^^ /(2p). At 77 ^ all baryons are accumulating 
near zero mass and a summation over all baryons gives us ~ instead of eq. JQ) ~ the following: 

$(r;, p) - $0 + ia< if + \bv^ + CTf'^p , (24) 

where C = "^icf / {2pi). The term p^^Pi is absorbed into $0 since it is independent on the chiral 
condensate 77. The typical momenta are of the order of the temperature, pi ^ T . Thus, Eq. (|24|l is valid 
in the region rj -^T^ . In the leading approximation the coefficient C can be considered as independent 
on the temperature at T ~ xjP^ . 

Due to the last term in Eq. (|24|l the thermodynamical potential always have a local minimum at 77 = 
since the condensate 77 is always non-negative. At small t < there also exists a local minimum at 77 > 0, 
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which is a solution of the equation: 



^ = (at + Brj'')7j+^Cp7j-^^^^0. (25) 




Figure 3: (a) Graphical representation of Eq. H25I) : "I" is the first term and "11" the second term (with 
the opposite sign) in the r.h.s. of the equation, (b) The thermodynamic potential l|24|l vs. the chiral 
condensate at a fixed baryon density p > 0. At low enough temperatures, T — Ti, the system resides 
in the chirally broken (hadron) phase. The first order phase transition to the quark phase takes place 
at Tph = T2 > Ti. At somewhat higher temperatures, T3 > Tph the system is in a mixed state. The 
temperature T4 = Tcont corresponds to a continuous transition to the pure quark phase, in which the 
thermodynamic potential has the form T5. 

At small enough baryon density p, Eq. (|25|l [visualized in Figure |2Ia)] has, in general, two roots, 
771 < 770 and r]2 > Jyo, where r/o = {—at/2tB)^/^ is the minimum of the first term in the right-hand side 
of Eq. 1)25(1 . The calculation of the second derivative d^^/drj^ shows that the second root 772 (if exists) 
corresponds to minimum of ^{rf) and, therefore, is a local minimum of $. The point 77 = 771 corresponds 
to a local maximum of the thermodynamical potential since at this point the second derivative is always 
non-positive. 

The thermodynamical potential $(7;,p) at (fixed) non-zero baryon density p has the form plotted in 
Figure |2^b). At low enough temperatures (curve Ti) the potential has a global minimum at 77 > and 
system resides in the chirally broken (hadron) phase. As temperature increases the minima at 77 = and 
at 77 = ?72 > becomes of equal height (curve T2 = Tph). At this point the first order phase transition to 
the quark phase takes place. At somewhat higher temperatures, T = T3 > Tph, the 77 > minimum of 
the potential still exist but $(r/ = 0) < $(772)- This is a mixed phase, in which the bubbles of the hadron 
phase may still exist. However, as temperature increases further, the second minimum disappears (curve 
T4 = Tcont)- This temperature corresponds to a continuous transition to the pure quark phase, in which 
the thermodynamic potential has the form T5. 

Let us calculate the temperature of the phase transition, Tph(/c), at non-zero baryon density p. The 
transition corresponds to the curve T2 in Figure |2Ib), which is defined by the equation <i>(?72,p) = $(77 = 
0, p), where ^2 is the second root of Eq. I(25|l as discussed above. The solution is 

/ V 3/5 / V 2/5 

1/2 

and the second minimum of the thermodynamic potential is at 772 — [4a (T*^"^ — Tph(p))/(5i3)] 

At a temperature slightly higher than Tph(/3) the potential is minimal at 77 = 0, but it has also an 
unstable minimum at some 77 > 0. The existence of metastable state is also a common feature of the first 
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order phase transition {e.g., the overheated Uquid in case of Uquid-gas system). With a further increase 
of the density p (at a given temperature) the intersection of the two curves in Figure 3(a) disappears and 
the two curves only touch one another at one point 77 = ^4. At this temperature a continuous transition 
(crossover) takes place. The corresponding potential has the characteristic form denoted as T4 in Figure 
3(a). The temperature T4 = T^ont is defined by the condition that the first (|25|) and the second derivatives 
of Eq. (|23I) vanish: 

/ \ 3/5 / V 2/5 

(I) , 

and the value of the chiral condensate, where the second local minimum of the potential disappears is 

/QN 1/2 

given by 774 — [2a{Tcont{p) — Tc ')/{5B)] . At temperatures T > Tcont(p) the potential has only one 
minimum and the matter is in the state with the restored chiral symmetry. Thus, in QCD with massless 
quarks the type of phase transition with the restoration of the chiral symmetry strongly depends on the 
value of baryonic density p. At a fixed temperature, T < Ti"', the phase transition happens at a certain 
critical density, pph- According to Eq. the critical density has a kind of a "universal" dependence 
on the temperature, Pph(r) oc [Ti^-* — T]^/"^, the power of which does not depend on the parameters of 
the thermodynamic potential, a and B. 
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Figure 4: The qualitative phase diagram at finite baryon density and temperature based on the analysis 
(a) without and (b) with indication of the approximate 2-nd order transition domain. 

The expected phase diagram is shown qualitatively in Figure 4(a). This diagram does not contain 
an end-point which was found in lattice simulations of the QCD with a finite chemical potential [25,26]. 
One may expect that this happens because in our approach a possible influence of the confinement on the 
order of the chiral restoration transition was ignored. Intuitively, it seems that at low baryon densities 
such influence is absent indeed: the deconfinement phenomenon refers to the large quark-anti-quark 
separations while the restoration of the of the chiral symmetry appears due to fluctuations of the gluonic 
fields in the vicinity of the quark. However, the confinement phenomenon dictates the value of the baryon 
size which can not be ignored at high baryon densities, when the baryons are overlapping. If the melting 
of the baryons happens in the hadron phase depicted in Figure 4(a), then at high enough density the 
nature of the transition could be changed. This may give rise in appearance of the end-point observed in 
Ref. [25,26]. 

The domain where the inequality \at\ ^ Cprf'/^, p 7^ is fulfilled, has specific features. In this 
domain the phase transition looks like a smeared second order phase transition: the specific heat has 
(approximately) a discontinuity at the phase transition point, ACp ~ a^Tc/2B. This statment follows 
from general theory [12]. At | at \^ Cpif/^ the last term in (17) may be neglected and we find for the 
entropy 

a$ „ IdA ^ 



In the phase above phase transition, rj = 0, S = Sq. Below phase transition 



(29) 



The specific heat Cp = T{dS/dT)p below the phase transition in the limit T — > Tc is equal 




(30) 



The correlation length increases as (T - rr) at T - rr 0. The latter arises if we include 
the derivative terms in the effective thermodynamical potential. The phase diagram with this domain 
indicated may look as it is shown in Figure 4(b). Note that the applicability of our considerations is 
limited to the region |T - T^"^\/T^"^ < 1 and low baryon densities. 

In the real QCD the massive heavy quarks (the quarks c, 6, t) do not influence on this conclusion, since 
their concentration in the vicinity of T w Ti°^ ~ 200 MeV is small. However, the strange quarks, the 
mass of which w 150 McV is just of order of expected Ti°\ may change the situation. This problem 
deserves further investigation. 

This work was supported in part by INTAS grant 2000-587, RFBR grants 03-02-16209. 
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